home *** CD-ROM | disk | FTP | other *** search
/ Inside Indy 1993 / Inside Indy 1993.iso / demos / CHEM / ribbon / data / 1hho_H.csh < prev    next >
Encoding:
Linux/UNIX/POSIX Shell Script  |  1993-06-23  |  625 b   |  26 lines
  1. #!/bin/csh
  2.  
  3. #  sample script to get heme subunits
  4. #   note the /bin/commands, as i have these aliased
  5.  
  6. # create termporary files with all FE and with no FE
  7. grep FE 1hho_H.pdb > fe.pdb
  8. grep -v FE 1hho_H.pdb > por.pdb
  9.  
  10. # loop over the numbers in 1hho  801..804
  11. foreach i ( 1 2 3 4 )
  12.     echo 80$i 80$i 11 > fe.in
  13.     echo 80$i 80$i 1 > por.in
  14.     pdb-range-sph fe.pdb 1hho_Fe$i.sph < fe.in
  15.     pdb-range-sph por.pdb tmp.sph < por.in
  16.     sph-bond tmp.sph 1hho_por$i.cyl
  17.     /bin/rm re.in por.in tmp.sph
  18.     end
  19.  
  20. #  pitch temporary files
  21. /bin/rm re.pdb por.pdb
  22.  
  23. # create ribbons files
  24. /bin/ls 1hho_Fe* > globin.atoms
  25. /bin/ls 1hho_por* > globin.bonds
  26.